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ABSTRACT 



Aims. We carried out a deep subarcsecond BRI imaging of the two middle-aged pulsars to establish their properties in the optical range 



o 



Methods. Astrometry and photometry methods are applied to identify the pulsars and to measure their fluxes. We also reanalyze archival 
»^ ESO/NTT and HST broadband data and find that some published fluxes for Geminga were estimated inaccurately. The resulting dereddened 
broadband spectra in the near-IR-UV range are analyzed and compared with available data from the radio through gamma-rays. 
Results. Both pulsars are detected at SlOcr level. Geminga is for the first time reliably detected in the / band with a magnitude of 25™ 10+0. 14. 
IL" The dereddened spectra of both pulsars are remarkably similar to each other and show significant flux increases towards the far-UV and 
. 5^ near-IR, and a wide flux excess in V-l bands. This suggests a multicomponent structure of the optical emission. The nonthermal power law 
component of the pulsar magnetospheric origin dominates in the most part of the optical range. For PSR B0656+14 it is compatible with a low 
5_j energy extension of the power law tail seen in hard X-rays. For Geminga the respective extension overshoots by a factor of 100 the nonthermal 
, optical flux, which has a less steep spectral slope than in X-rays. This implies a spectral break at a photon energy ~ 1 keV. The flux increases 
towards the far-UV are compatible with contributions of the Rayleigh- Jeans parts of the blackbody components from whole surfaces of the 
neutron stars dominating in soft X-rays. The V-I excess, which is most significant for PSR B0656+14, suggests a third spectral component of 
still unidentified origin. Faint, a few arcseconds in size nebulae extended perpendicular to the proper motion directions of the pulsars, are seen 
around both objects in our deepest / band images. They can be optical counterparts of the bow-shock head of Geminga and of the tentative 
pulsar wind nebula of PSR B0656+14 observed in X-rays. 
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1. Introduction 

The two nearby middle-aged pulsars, PSR B0656-fl4 and 
Geminga (PSR J0633-(-1746), belong to a small group of ro- 
tation powered isolated neutron stars (NSs) whose emission is 
detected in a wide spectral range spanning from radio through 
gamma-rays. Their parameters, such as the period P and its 
derivative P, characteristic age t, spin-down luminosity E, 
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* Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 



magnetic field B, are very similar (Table This is reflected 
in similarities of observational manifestations of both pulsars. 

Multiwavelength spectra of PSR B0656-I-14 and Geminga 
consist of two main components. The first one is a nonther- 
mal component which is believed to be generated in magne- 
tospheres of the NSs. It dominates the pulsar emission in al- 
most the whole observed range and its spectrum is described 
by a power law, albeit with different spectral indices in differ- 
ent spectral domains. Multiwavelength properties of this com- 
ponent are important for the study of not yet clearly understood 
emission mechanisms in magnetospheres of NSs. The second 
component dominates in soft X-ray and EUV ranges where it 
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is seen as a strong excess over the nonthermal background. 
It is well described by a blackbody spectrum and this excess 
is thought to be associated with thermal emission form cool- 
ing surfaces of the NSs (Becker & Triimper 1997 Zavlin & 
Pavlov 1200431 De Luca et al. l2005> . The study of this compo- 
nent is of a crucial importance to evaluate the NS surface tem- 
perature, to compare it with NS cooling theories and to obtain 
information on still poorly known properties and equation of 
state of super-dense matter in interiors of NSs (e.g., Yakovlev 
& Pethick 2004 1. Inclusion of an additional high temperature 
thermal component assumed to be originated from hot pul- 
sar polar caps appears to improve significantly spectral fits in 
X-rays (Koptsevich et al. 2001, Pavlov et al. 2002; Zavlin & 
Pavlov l2()()4al De Luca et al. l2005l . 

There are, however, noticeable differences between the two 
objects. Geminga is known as one of the brightest gamma- 
ray sources on the sky (Bignami & Caraveo fl996> . while PSR 
B0656-I-14 is fainter and only marginally detected in gamma- 
rays (Ramanamurthy et al. I1996t . In the radio range PSR 
B0656-I-14 is a stable source with properties to be typical for 
an ordinary radio pulsar In contrast, Geminga has been consid- 
ered for a long time as a "radio-quiet" NS, until it was detected 
at rather low (only ^ 102.5 MHz) frequencies (Kuz'min & 
LosovskiiTggT", Malofeev & MalovTO^T, Shitov & Pugachev 
[1998 1. Geminga has unusually steep radio spectrum (Burderi 
et al. 1999), and unlike most of radio pulsars its pulse profile is 
very unstable (Kassim & Lazio il999> . 

PSR B0656+14 has a twice lower transverse velocity, 
Vi ^ 60 km s ' (Mignani et al. 2000; Brisken et al. 2003 1, 
than Geminga, v^ ^ 122 km s ' (Caraveo et al. ,1996 1, and 
both pulsars move in different environment, as indicated by 
their dispersion measure values (cf. Table 1). PSR B0656H-14 
sits near the center of the Monogem ring (Nousek et al. 1981 ), 
that is a bright ~ 10^ yr expanding supernova remnant (SNR) 
likely produced at the same supernova explosion as the pul- 
sar (Thorsett et al. i2003il . Geminga has not been associated 
with any known SNR, although based on its proper motion 
and age a birth place was proposed to be located inside the 
Cas-Tau OB or the Ori OB la stellar associations (Pellizza et 
al.|2005 1. Recent XMM-Newton observations of this pulsar in 
X-rays revealed a faint bow-shock structure whose tails extend 
up to 2' from the pulsar and are elongated along direction of 
its proper motion (Caraveo et al. 2003 1. The structure has not 
yet been identified in the optical range. No reliable signs of 
a bow-shock nebula have been detected around more distant 
and slowly moving PSR B0656-I-14. Some hints of a compact, 
3"-4", faint pulsar wind nebula (PWN) surrounding the pul- 
sar were reported based on the X-ray images obtained with the 
Chandra/LETG/HRC-S (Marshall & Schulz l2002l and with the 
Chandra/ ACIS-S3 (Pavlov et al. 2002 1, but they have not yet 
been confirmed by deeper exposures and have not been identi- 
fied in other spectral domains. 

Both pulsars are firmly identified in the optical range 
by an accurate positional coincidence of the optical counter- 
parts with the pulsars, by comparison of the proper motion 



Table 1. Parameters of Geminga and PSR B0656-H14 (from 
Taylor et al. ll9931 unless specified otherwise). 



Geminga 



B0656+14 



Observed 

a (J2000)° .... 6*33'"54n530(28)* 

6 (J2000)° .... HMe' 12'.'909(40)'' 

(masyr-'). 138(4/ 

(masyr-') . 97(4)'' 

;r°(mas) 6.36(1.74)'' 

f 195° 134 

b' 4° 265 

P (ms) 237 

P(10-'^) 10.97 

DMf {cm-^ pc) 2.9« 

Derived 

cPipc) l57CJ,y 

T* (kyr) 342 

logS(G) .... 12.21 

log£:(ergs"') 34.57 



6*59'"48n472(7)'' 

14''14'21'.'160(10)'' 

45.47(65)" 

-2.40(29)'' 

3.47(36)" 

201° 108 

8° 25 8 

385 

55.01 

14.02 



288(^33)^ 

111 

12.67 

34.58 



in radio and optical ranges (Caraveo et al. 119961 Mignani et 
al. EOOO. Brisken et al. 2003 1, and by detection of pulsations 
with the pulsar period in the optical (PSR B0656-I-14; Shearer 



" coordinates, proper motion, parallax, and parallax-based 
distance; hereafter the numerals in parentheses are Icr 
uncertainties referring to last significant digits quoted; 
Caraveo et al. 1998 (epoch 1995.21); 

" Brisken et al. 2003j 

'' Cai-aveo et al. 1995] 

'' galactic coordinates; 

^ dispersion measure; 

11 Malofeev & Malov lTW/l 
characteristic age P/2P. 



et al. nWTl Kern et al. '5005' Geminga: Shearer et al. 

and in UV (PSR B0656+14: Shibanov et al. l2005l Geminga: 

Kargaltsev et al. 2005 1. 

Due to faintness of the pulsar counterparts the spectral 
information in the optical is restricted mainly by the broad- 
band photometry. Published spectral observations of Geminga 
(Martin et al. .1998.) are too noisy to compete with the broad- 
band data. The broadband far-UV (FUV) spectra of both 
pulsars appear to be dominated by the Rayleigh-Jeans (RJ) 
part of the blackbody spectrum extrapolated from soft X-rays 
(Mignani et al. 1998 Koptsevich et al. 2001 hereafter WMT\ . 
For Geminga this is supported also by FUV spectral data re- 
cently reported by Kargaltsev et al. (2005 1. At longer wave- 
lengths there is a significant flux excess over the RJ extrapola- 
tion in the spectra of both pulsars. 

For each of the pulsars these excesses have been interpreted 
in a different way. For PSR B0656H-14 the spectral energy dis- 
tribution from the optical to near-IR (NIR) ranges was shown 
to be compatible with the low energy extension of a high en- 
ergy nonthermal spectral tail detected in X-rays (K|2001 1. It 
suggested possible common origin of the nonthermal spectral 
component in the optical and X-rays. Early attempts of detec- 
tion of Geminga in NIR were unsuccessful and the measured 
upper limit in the / band suggested a strong spectral dip in this 
band followed by a prominent excess in the V band (Mignani et 
al. [1998). The shape of the excess reminded a broad emission 
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Table 2. Log of observations of Geminga and PSR B0656+14 
with the Subaru. 



PSR 


UT 


Exp. 


Band Seeing" sec z 




21 Jan 2001 


sec 


arcsec 



Geminga 05:17:26.9 332.3 / 0.66 1.56 

Geminga 05:30:12.6 300.0 / 0.66 1.47 

Geminga 05:42:29.2 300.0 / 0.70 1.39 

Geminga 05:51:53.6 300.0 / 0.65 1.35 

Geminga 06:01:30.5 300.0 / 0.68 1.30 

Geminga 06:31:27.7 300.0 / 0.55 1.19 

Geminga 06:40:32.8 300.0 / 0.55 1.16 

Geminga 06:49:38.0 400.0 / 0.55 1.13 

Geminga 07:00:22.8 400.0 / 0.60 1.11 

Geminga 07:11:07.8 400.0 / 0.58 1.09 

Geminga 07:21:52.8 400.0 / 0.60 1.07 

Geminga 07:32:38.2 400.0 / 0.60 1.05 

Geminga 07:43:23.9 400.0 / 0.60 1.04 

Geminga 08:07:55.3 400.0 / 0.61 1.01 

Geminga 08:26:36.9 400.0 / 0.62 1.00 

B0656+14 08:46:10.7 400.0 / 0.65 1.01 

B0656+14 08:56:55.8 400.0 / 0.65 1.01 

B0656+14 09:07:40.9 400.0 / 0.73 1.00 

B0656+14 09:18:26.0 400.0 / 0.67 1.00 

B0656+14 09:29:10.9 400.0 / 0.65 1.01 

B0656+14 09:39:55.9 400.0 / 0.63 1.01 

B0656+14 10:22:05.8 510.0 R 0.60 1.04 

B0656+14 10:34:41.0 510.0 R 0.60 1.06 

Geminga 10:47:15.8 510.0 R 0.61 1.13 

Geminga 10:59:58.6 510.0 R 0.68 1.17 

Geminga 11:12:33.8 510.0 R 0.61 1.21 

Geminga 11:25:09.0 510.0 R 0.60 1.25 

B0656+14 11:54:36.9 600.0 B 0.65 1.28 

B0656+14 12:08:42.3 600.0 B 0.66 1.35 

B0656+14 12:22:47.8 600.0 B 0.69 1.43 
° FWHM of the stellar profile in the individual images. 

feature over the RJ continuum that has been interpreted as an 
ion cyclotron emission from a hot non-uniformly magnetized 
plasma covering the pulsar polar caps (Jacchia et al. 1999 1. 
However, later Geminga was detected in the NIR (Harlow et 
al. 0[998). Preliminary data analysis ({ 3200 1> revealed a flux 
growth towards longer wavelengths, as it is in case of PSR 
B0656+14. This does not quite agree with the cyclotron inter- 
pretation. In addition, Kurt et al. (2001 1 reported on a tentative 
detection of Geminga in the / band at a flux level, which is by a 
half order of magnitude higher than the upper limit known pre- 
viously. The reported flux value is compatible with a smooth 
connection between the V and NIR bands, that has made the 
cyclotron interpretation even more questionable. 

In any case, published broadband optical data are still too 
uncertain to reject or confirm surely the suggested interpreta- 
tions and to understand whether these two middle-aged pulsars 
have similar or diff'erent properties in the optical range. The 
uncertainties do not allow us to study also significance of the 
spectral difference between the middle-aged and other pulsars 
of diff'erent ages detected in optical range ( H20011 Zharikov et 
al.p002, 2003, Shibanov et al. 2003 1. The main source of the 
flux errors in case of PSR B0656-I-14 is a contamination of the 
pulsar flux by a background extended object in ^ 1'.' 1 from the 



pulsar (K I2001t . For about twice fainter Geminga, there are no 
contaminating objects but deeper exposures are necessary to 
establish its optical spectral energy distribution with a higher 
accuracy. 

In this paper we report on a subarcsecond deep broadband 
imaging of the Geminga and PSR B0656-H14 fields with the 
Subaru telescope. Preliminary results have been reported by 
Komarova et al. ( 2003 ). Here we analyze our data together with 
the available NIR-optical-UV data from the HST, ESO/NTT, 
and 6m BTA telescope, as well as with available multiwave- 
length data from the radio through gamma-rays. In Sect. |2 we 
describe the observations, data reduction, astrometric and pho- 
tometric referencing. The results are presented in Sect. |3l and 
discussed in Sect.|3 

2. Observations and data analysis 

2.1. Subaru observations and data reduction 

Two fields containing Geminga and PSR B0656-H14 were 
observed on January 21 2001 with the wide field camera 
Suprime-Cam at the primary focus of the Subaru telescope. 
The Suprime-Cam (Miyazaki et al. 1998i was equipped with 
ten MIT/LL 2048 x 4096 CCDs arranged in a 5 x 2 pattern 
to provide 34' x 27' FOV with the pixel size on the sky of 
0.201" X 0.201". The pulsars were exposed in one of the CCD 
chips, si(506s, with the Gain=1.17, making use of the B, R, and 
/ filters with throughput to be close to the Johnson-Cousins sys- 
tem. Sets of 5-10 min individual dithered exposures were ob- 
tained in each of the filters to get rid of cosmic rays at further 
data analysis. Several individual exposures were removed from 
the consideration due to problems with the telescope guiding 
system. Thus the total exposure time used for the subsequent 
analysis of Geminga was 2040 s in the R and 5332.3 s in the 
/ bands. For a brighter PSR B0656+14 it was 1800 s, 1020 s, 
and 2400 s in the B, R, and / bands, respectively. The observa- 
tional conditions were rather stable during the night with seeing 
varying from ~ 0'.'6 to ~ 0'.'7. A Log of observations is given 
in Table 13 The PG 1047 +003 standard field of Landolt (.19921 
used for the photometric calibration was observed at the end of 
the night. 

Standard data reduction including bias subtraction and flat- 
fielding was performed making use of the MIDAS and IRAF 
software. In order to get rid of the remaining cosmic rays and 
to decrease possible contamination of PSR B0656+14 by a 
nearby extended object we combined the individual dither im- 
ages of the pulsar field in the B and / bands making use of the 
ditherll package by Fruchter et al. 1 19981 

2.2. Astrometric referencing of the Subaru images 

The positions and proper motions of Geminga and PSR 
B0656-I-14 are well established (see Table^J. To determine the 
expected coordinates of Geminga at the epoch of the Subaru 
observations (2001.054) its position at the epoch of the ob- 
servations with HST/WFPC2 (1995.21) (Caraveo et al. [T996l 
was used as a reference point and then corrected for the proper 
motion. The same was done for PSR B0656+14 using a refer- 
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Table 3. Positions of Geminga and PSR B0656+14 at the epoch of the Subaru observations, 2001.054. 



Object 


Expected 


Observed 




Q'2000 "52000 


Q'2000 '^2000 


Geminga 


06"33"'54.2067(32) 17° 46' 13'.' 475(46) 


06"33"'54'.219(15) 17°46' 13'.' 56(22) 


PSRB0656+14 


06'' 59"'48n503(7) 14° 14'2 1'.' 1 575( 10) 


06''59"'48M52(14) 14°14'21'.'41(21) 




48.8 48.4 6:59:48.0 47.6 48.8 48.4 6:59:48.0 47.6 48.8 48.4 6:59:48.0 47.6 

Fig. 1. 20" X 20" fragments of the Subaru images of the PSR B0656+ 14 field in the B, R, and / bands. Units of RA (horizontal 
axis) are hours, minutes and seconds, and units of Dec (vertical one) are degrees, arcminutes and arcseconds. A box in the right 
panel marks the region which is magnified in Fig.|6land discussed in Sect. 13. 21 



ence point from Brisken et al. Il2003> . based on the radio obser- 
vations with the VLBA. The expected pulsar coordinates are 
given in Table |3] with errors accounting for the proper motion 
and reference point uncertainties. 

The astrometric referencing of the Subaru images was done 
making use of the USNO-Bl.O catalog (Monet et al.TOOZi and 
IRAF tasks ccmap/cctran. For the astrometric transformation 
of the combined images of the Geminga field we used the po- 
sitions of six reference stars'. Formal nns errors of the astro- 
metric fit for the RA and Dec were ^ U.'Ql and ^ (^.'06, respec- 
tively, for both the R and / images. Seven reference stars^ were 
used in case of PSR B0656+14. For each of the B, R, and / im- 
ages formal nns errors were (K'OS and a; (X'OS, for the RA and 
Dec, respectively, which is less than the nominal USNO-Bl.O 
catalogue accuracy of (J.'IQ. Combining the rms errors and the 
catalog uncertainty we estimate the total 1 cr accuracy of our as- 
trometric referencing of all the images as 0'.'21 in both RA and 
Dec. The optical counterparts to both pulsars are clearly seen at 
the expected positions (see Fig.^and Fig.O. The expected and 
observed positions (Table O coincide with each other within 
the errors accounting for the centroid and astrometric referenc- 
ing uncertainties. 



' USNO-Bl.O stars used for the astrometric transformation of the 
Geminga field: 1077-0137329, 1078-0139943, 1077-0137377, 1078- 
0139792, 1078-0139780 and 1077-0137166 

^ USNO-Bl.O stars used for the transformation of the PSR 
B0656+14 field: 1042-0123510, 1042-0123514, 1042-0123581, 
1042-0123600, 1042-0123481, 1042-0123424 and 1042-0123352 



2.3. Photometric referencing of the Subaru images 

Photometric calibration was carried out using three standard 
stars from the Landolt field PG1047-H003 (Landolt [T992t . ob- 
served the same night and five secondary standard stars from 
the PSR B0656+14 field (Kurt et al. 1998 1. The atmospheric 
extinction factors in the B, R, and / bands were derived from 
the variation of the count rates of field stars with the air- 
mass changing during our observations (see Table|2j. The val- 
ues obtained kg = 0™18 + 0.02, ku = 0^117 + 0.021, k, = 
0^042+0.020 are consistent with the mean values at the Subaru 
site (Y. Komiyama, private communication). Insignificant de- 
crease of the extinction, within Icr level, was noticed from the 
beginning to the end of observations. Based that the following 
transition equations related to the observational night were ob- 
tained: 

B-b^ 28.283(20) + 0.094(10) -(b-r) 
R-r = 28.122(29) - 0.023(14) -(b-r) 
R-r^ 28.070(44) - 0.138(47) ■ (r - /) 
I-i^ 27.340(38) - 0.152(55) ■ (r - i) 

where b, r, i are instrumental magnitudes and B, R, I are mag- 
nitudes in the Johnson-Cousins photometric system. The eiTors 
of the transformation coefficients shown in parentheses include 
the extinction factor and standard instrumental magnitude un- 
certainties. 

2.4. The data from the ESO/NTT and HST archives 

In our analysis we used also the I-band images of the Geminga 
field obtained with the ESO/NTT (Bignami et al. ,1996j . kindly 
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Object 


Band 


seeing" 


5/^ 


d" 


Source"^ 


kl cosCZ/ 




Mag^' 


Flux 






pix. 




pix 


-2.5 Log(counts/s) 


mag 


mag 






PSR B0656+14 


B 


4.5 


22 


4 


-1.68 


0.25 


1.12 


25.27(7) 


0.313(20) 




R 


3.4 


18 


4 


-2.50 


0.07 


0.89 


24.65(8) 


0.417(30) 




I 


4.1 


11 


4 


-1.90 


0.05 


0.95 


24.52(13) 


0.371(42) 


Geminga 


R 


3.7 


15 


4 


-1.63 


0.14 


0.84 


25.53(11) 


0.185(17) 




I 


3.9 


10 


4 


-1.37 


0.05 


0.87 


25.10(14) 


0.217(26) 



° FWHM of the stellar profile in the combined images. Corrections for the atmospheric extinction. 

* Diameter of the optimal photometrical aperture. '' Correction factors for the PSF. 

Instrumental magnitudes, m = Source - kj cos(Z) - Sfmd. ^ Johnson-Cousins magnitudes. 




17:46:05 

-k 



54.fi 54.4 6:33:54.0 53.6 

• Subaru, 21 jaif''j"ob^-'|:'^^^ 
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Fig. 2. 20" X 20" fragments of the Subaru images of the 
Geminga field in the R {top) and / (bottom) bands. Axes ori- 
entation and label notations are the same as in Fig. 1. 



provided by R. Mignani. Five 20 min science exposures were 
obtained on January 31 1995 using the SUSI^. The reduction 
of these data was rather complicated due to the absence of a 
standard "flatfield" image in the archive. To produce an arti- 
ficial flatfield we subtracted stars from each science exposure 
and combined these images into a median one. It was slightly 
improved in the regions around the star positions making use of 
the MODIFY/AREA task of the ESO/MIDAS package, smoothed 
with the 5x5 pixel average and used as the flatfield in a standard 
data reduction. The problems with the flatfield and the fact that 

^ Observational conditions and other details can be found in 
Bignami et al. 1996 



the individual science images were taken without shifts enough 
to remove accurately all artifacts led to the appearance of arti- 
ficial extended objects around the position of Geminga in the 
final image (Fig.|3] left panel), precluding from a reliable iden- 
tification of the pulsar (see Sect. 3.1. for details). 

The NIR-optical-UV data from the HST archive obtained 
wifli the HST/NICMOS in the FllQW and F160W bands 
(HST analogs of J and H\ Harlow et al. TWl and with 
the HST/WFPC2 in the F555W band for Geminga, with the 
HST/FOC in flie F\95W , F342W, F430W_bands for Geminga 
(Bignami et al. [T996I Mignani et al. and for PSR 

B0656-H 14 (Pavlov et al. [T996inW7t were also reanalysed. The 
data reduction and analysis of the NICMOS images obtained at 
the epoch of 1998.22 were performed as described in K2001P. 
The photometry of the WFPC2 images processed by pipeline 
was performed using the algorithm described in the WFPC2 
Data analysis Tutorial''. The reduction of the FOC images was 
done in the same manner as described in Pavlov et al. (19961 
and ll997t making use of the FOC Instrument Handbook (v6)^. 



3. Results 

The optical counterparts of both pulsars are firmly detected in 
all Subaru images with a signal to noise ratio S /N > 10. The in- 
strumental magnitudes of the counterparts were measured for a 
range of circular aperture radii of (1 - 3) CCD pixels centered at 
the source positions. They were corrected for the point spread 
function (PSF) of bright stars. The background was measured 
over the annulus with the inner and outer radii r,„ = 5 and 
four -10 pixels. Within the measurement errors the results for 
different apertures coincided and a 2 pixel radius, or 4 pixels 
in the diameter that is comparable with the seeing of about 3-4 
pixels (cf. Table 2), was adopted as an optimal to get a max- 
imum S/N. Using Eq (1) Johnson-Cousins magnitudes were 
derived and transformed into the fluxes in physical units using 
zero-points provided by Fukugita et al. 119951 The results are 
given in Table |3 The errors include statistical errors of the in- 



archive.stsci.edu 

5 See also the HST Data Handbook for NICMOS, 
stsci.edu/hst/nicmos/documents/handbooks/ 

^ Baggett et al. 2002, in HST WFPC2 Data Handbook, v4, ed. B. 
Mobasher, Baltimore, STScI 

' _stsci.edu/ftp/instrument news/FOC/Foc handbook| 
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Fig. 3. Left: A fragment of the Geminga field image obtained with the ESO/NTT in the / band (Bignami et al. I1996> . The 
contours of the Subaru / band image are overlaid and the dashed-line arrow shows the direction of the pulsar motion; the cross 
marks the position of the pulsar at the epoch of the ESO/NTT observations. Middle and right: Near-infrared images of the same 
field obtained with the HST/NICMOS in the Fl IQW and F16QW bands, respectively (Harlow et al. ll998> . Boxes mark regions 
with Geminga, they are magnified in Fig.|3 



strumental magnitude measurements, the errors of photometric 
referencing, and possible atmospheric extinction variations. 

3.1. Geminga 

3.1 .1 . Identification in the I band 

The object of /a;25™l detected with the Subaru at the expected 
Geminga position is by a magnitude brighter than the upper 
limit />26™4 reported by Mignani et al. ( 1998 1 based on previ- 
ous / band observations with the ESO/NTT. Such a deep limit 
is puzzling since it is comparable with the 3cr detection limit of 
our Subaru observations, I> 26'."8, obtained with similar expo- 
sure time (1.5 h vi 1.4h of NTT), but with a much higher tele- 
scope collective area (8 m v,s- 3.6 m of NTT) and at better seeing 
conditions (0'.'6-a.'8 vs O'.'9-l'.'l of NTT). Using the NTT/SUSI 
Exposure Time Calculator (ETC)^ one can easily estimate that 
at the respective conditions the NTT image should be by a mag- 
nitude less deep than the reported limit. This is in agreement 
with an upper limit of />25™4, which is presented in the text of 
the paper by Bignami et al. f 1996^ describing the NTT obser- 
vations of Geminga. However, this reasonable limit was likely 
incorrectly reproduced in Figs. 3 and 4 of the same paper, and 



later on in Mignani et al. J1998I I. 

To solve this puzzle we inspected the archival NTT im- 
age, which is shown in the left panel of Fig.|3l The expected 
position of Geminga at the epoch of the NTT observations, 
1995.084, is marked by a cross and the proper motion of the 
pulsar is shown by a dashed arrow. Several background-like ob- 
jects are seen in this image in a few arcseconds from the posi- 
tion of Geminga, while the pulsar itself can be hardly identified 
in this image. A thorough comparison with the Subaru /-band 
image {bottom panel of Fig. |2}, whose contours are overlaid, 
shows that these objects are likely to be artifacts caused by a 
poor data reduction in the absence of the standard Hatfield im- 



age in the archive (cf. Sect. 2.4). This is further strengthened by 
the fact that there are no background objects in this region de- 
tected neither in our Subaru R band {top panel of Fig.|2} nor in 
the HST F\ lOW band images (Fig.|3]E} which are adjacent to 
/ band in the spectral energy distribution and where Geminga is 
clearly identified. Our rough estimate of the 3cr detection limit 
of the NTT image shown in Fig. Ogives /> 24™4. This is by a 
magnitude less deep than the above ETC and those of Bignami 
et ai. J1996> estimates and can be caused by the data reduc- 
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Fig. 4. 4" X 4" images of Geminga in the FllQW {left) and 
F160W {right) bands. The circle annulus in the left panel mark 
the expected Geminga positions at the epochs of I-band obser- 
vations with the NTT (1) and Subaru (2), respectively, while the 
long arrow shows the pulsar path calculated based on the proper 
motion measurements of Caraveo et al. (J 996 ). Contours of the 
Subaru image in / band are overlaid in the right panel where the 
arrow shows the shift of Geminga due to its proper motion over 
3.2 yr between the HST and Subaru observations; and Icr- 
ellipse mark the expected positions of the pulsar at the epochs 
of the HST and Subaru observations, respectively. 



tion problems mentioned above. However, it is hardly possi- 
ble to get a two magnitude deeper flux limit, as in Mignani 
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et al. <1998> . by any improvement of the NTT data reduction. 
Therefore, it has been overestimated at least by a magnitude for 
unknown reasons. On the other hand, we note that the object 
detected at the Geminga position in the Subaru / image could 
be likely detected also in the NTT image with S /N ~ 4, pro- 
vided a straightforward NTT data reduction is possible. This 
would be consistent with the significance of a marginal detec- 
tion of Geminga in the / band with the 6m BTA telescope re- 
ported by Kurt et al. ( 1200 1> . 

To insure further our / band and NIR identification of 
Geminga we estimated the shift between the counterpart po- 
sitions in the HST/NICMOS and Subaru images (see Fig.|4j. It 
was found to be 0'.'45±0'.'22. Despite the only 2cr significance 
of this estimate it is in a good agreement with the shift direction 
(marked by arrows in Fig.|3 and with its value of O'.'31+O'.'Ol, 
which are expected from the Geminga proper motion (Caraveo 
et al. 1 19961 Tabled at the time base of about 3.2 yr separating 
the HST and Subaru observations. We have not found also any 
reliable faint red background source on the path of the pulsar 
proper motion, which would prevent from correct identifica- 
tion of Geminga in the / band at epochs after the NTT observa- 
tions, as was warned by R. Mignani (private communication). 
This allows us to make a conclusion on a firm identification of 
Geminga in the /, as well as in R, and in the two NIR bands. 

3.1 .2. Broadband fluxes from the NIR through FUV 

The results of our I-R and NIR photometry of Geminga are pre- 
sented in Table |5l where we have collected also most accurate 



Table 5. Geminga in the NIR-optical-NUV-FUV range. 



Filter/ 
Telescope 


Log V 
Hz 


observed 
//Jy 


dereddened" 
//Jy 


NIR 


F1601V/HST 
FUOW/HST 


14.273(+^^) 
14.436(!|38) 


0.230(25)''-" 
0.166(21)*-^ 


0.233(25) 
0.170(21) 


optical 


//Subaru 
i?/Subaru 
F5551V/HST 
F5551V/HST 


14.571(!^^) 
14.658(+55) 

14.734(+g) 
14.734(+g) 


0.172!|^^ 
0.170(20)'' 
0.146(13)^ 


0.20 

0.181![^ 
0.182(21) 
0.156(14) 


F430M'/HST 


14.882(:«) 


0.146(11)'- 


0.160(12) 


NUV 


F3421V/HST 

F25SRF/HST 

F195W/HST 


14.944(!*0 
15.116(!jf) 
15.104(!^5'') 


0.201(10)^ 
0.261(26)'' 
0.288(30)" 


0.224(11) 
0.315(31) 
0.344(36) 


FUV 


G140L/HST 


15.322(84) 


0.455(29)'' 


0.546(35) 



" E(B-V)=0.023^ ' this work; 

* Harlow et al. [l998l '' Kargaltsev et al. l2005l using ACS 
and STIS/MAMA NUV and FUV. 



and partially updated results of available pulsar observations in 
other broad bands. The flux in the R band is in a good agree- 
ment with what was published by Mignani et al. ( .1998 ) , but 



the accuracy of the Subaru observations is considerably higher. 
The result in the / band confirms a tentative detection reported 
by Kurt et al. ( 2001 ), but at much higher significance level. The 
NIR fluxes are in agreement with the preliminary results pub- 
lished by Koptsevich et al. ( 1200 U and Komarova et al. (12003> . 

Our reanalysis of the four HST/WFPC2 datasets obtained 
at different epochs in the F555W band in the period of 1994- 
1 995 yield the Geminga fluxes in the range of 0. 1 36-0. 161 fiiy 
depending on the set. The flux values in separate sets are con- 
sistent with each other and do not show any variation of the 
Geminga emission with time above 3cr level. In Table |5l we 
show a mean F555W flux over these sets. It is also consis- 
tent with the value obtained recently with the HST/ ACS in the 
same band (Kargaltsev et al. 2005, cf. Table |5jl but a factor of 
1.7 smaller than early results (Bignami et al. 1996, Mignani et 
al. 1998 1. This provides independent confirmations of a signifi- 
cant overestimation of the Geminga flux in the F555 W band in 
the previous publications, as it has been noticed by Kargaltsev 
et al. ( 2005 ). At the same time, our reanalysis of the HST/FOC 
archival data obtained in F430W, F342W, and F195W bands 
shows a less significant, about 2cr (or 10%-20% of the flux 
value), but systematic underestimation of the previously pub- 
lished fluxes. In Table|5lwe present the updated flux values in 
these bands. 



3.1 .3. Possible bow-shock structure 

Observations of the Geminga field in X-rays with the XMM 
have revealed a faint bow-shock structure around this pulsar 
whose tails are extended up to 2. '5 behind the moving NS 
(Caraveo et al. 2003). Bow-shock structures produced by in- 
teraction of supersonically moving NSs with ambient medium 
have been detected around several pulsars presumably in the 
H(j, emission (e.g., the Guitar nebula around PSR B2224-I-65; 
Cordes et al. ll993> . where they are usually easier detected than 
in X-rays (e.g., Romani et al. 1997). However, no signs of any 
Ha counterpart of the Geminga shock structure were detected 
at deep imaging carried out in naiTow //„ bands with the VLT 
(Caraveo et al. l2003t and 6m BTA (Komai'ova et al. l2005t . 
A faint Ha counterpart, if exists, is likely hidden in a rather 
strong extended background emission in the Ha line from a low 
Galactic latitude region where Geminga is located (b~4'?7). 

At the same time, as in case of the Crab pulsar nebula 
(Hester et al. 2002), the brightness of the shocked structures 
of Geminga could be higher in the optical continuum instead 
of lines. In this case they have to be more easier detected with 
broadband observations. A careful inspection of the Subaru im- 
ages has allowed us to find possible indications of the optical 
counterpart of the X-ray structure in the deepest / band image. 
A faint arc-like feature resembling a head of the bow-shock 
can be seen in Fig.|5l Its contours are shown in the right panel 
of Fig. 0] The arc length is ~ 1'.'9 and it extends perpendic- 
ular to the pulsar path. There are no signs of this structure in 
our less deep R band image as well as in the archival HST 
images. Further deep imaging of the field is needed to under- 
stand whether the arc-like structure is indeed a counterpart of 
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Fig. 5. A magnified fragment of the Subaru / band image of 
Geminga demonstrating an arc-like structure resembling a head 
of a bow-shock due to the pulsar supersonic motion in the in- 
terstellar matter (cf. contours in Fig|^. The motion direction is 
shown by arrow, the pulsar position is marked by a cross. 

the bow-shock seen in X-rays, a faint background extended ob- 
ject, or an artifact. 

3.2. PSRB0656+14 

3.2.1 . Broadband fluxes from the NIR through NUV 



tion and to measure for the first time the fluxes in these bands 
with the accuracy of better than 10%. The results of our pho- 
tometry of the pulsar are shown in Table |6] where the most 
accurate and informative data in other broad bands are also col- 
lected^. 




Fig. 6. Magnified fragment of the Subaru / band image of PSR 
B0656+14 demonstrating possible existence of a PWN or a 
bow-shock structure around the pulsar (extended object in the 
center of the image). Contours are overlaid to underline the 
shape of the structure; the arrow shows the direction of the pul- 
sar proper motion, the pulsar position is marked by a cross 



Table 6. PSR B0656+14 in the NIR-optical-NUV range 



Filter/ 


Log V Fy 


Fy 


Telescope 


observed 


dereddened" 




Hz fi]y 




NIR 


F1871V/HST 


14.204^30) 0.779(97)* 


0.789(98) 


FieOW/HST 


14.273(;58) 0.575(41)* 


0.584(42) 


FUOW/HST 


14.436C|^^) 0.336(30)* 


0.346(31) 


optical 


//Subaru 


14.511(1^) 0.31l(l^,y 


0.390(+^^) 


iJ/Subaru 


14.65SC_11) 0.417(+«)- 


0.447(+«) 


F5551V/HST 


14.734(+g) 0.393(24)* 


0.428(26) 


B/Subaru 


14.829(!J2) 0.3l3Cl^y 


0.350(!^^) 


F4301V/HST 


14.882(+^5) 0.264(29)'' 


0.300(33) 


NUV 


F3421V/HST 


14.944(!^^) 0.307(27)" 


0.353(31) 


GRISM/nST 


15.071(+]™) 0.354(31)" 


0.429(38) 


F195W/HST 


15.104(+jf) 0.349(43)'' 


0.438(54) 


° E(B-V)=0.03; 


'' Pavlov etal. 1997 




* kEooTI 

this work; 


" Shibanov et al.l2005l 





Previous photometric measurements of the pulsar fluxes in 
the B, R, and, especiaUy in the / band (Kurt et al .l998. K.2001J 
were rather uncertain due to poor seeing conditions and a con- 
tamination from a faint extended background object o2 which 
sits only in I'/l from the pulsar (cf. Fig.Q. Perfect seeing dur- 
ing the Subaru observations allowed us to avoid the contamina- 



3.2.2. Possible PWN structure 

As in case of Geminga we find also a hint of a faint com- 
pact nebula around PSR B0656-i-14 possibly associated with 
the pulsar activity. This is seen in the deepest image obtained in 
the / band (Fig.|6|l. The detected structure is elongated approx- 
imately perpendicularly to the direction of the pulsar proper 
motion. Its maximal extension is ~ 175. If this is not a back- 
ground object or an artifact, it is the second indication of the 
presence of a faint compact PWN or a bow shock nebula around 
PSR B0656+14 in the optical range. The first one was men- 
tioned by Koptsevich et al. (2001i who found a knot-like struc- 
ture seen at a 3cr level above the wings of the pulsar spa- 
tial profile at ~ 0'/04 from the pulsar position in the archival 
HST/FOC/F343W image. Knot-like structures are typical at- 
tributes of the Crab PWN. However, in case of PSR B0656-t-14 
both indications still have to be considered with caution and 
need to be confirmed by deeper observations since we found 
no signs of the nebula in any other available optical images 
of the pulsar. However, all these images are less deep than the 
Subaru / band image. 



' As in case of Geminga (Sect. 3.1), using archival data we have 
examined the published fluxes of PSR B0656+14 obtained with dif- 
ferent HST instruments. No significant discrepancies from the results 
as they are listed in t EOOTI were found. 
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4. Discussion 

Reliable detection of Geminga and PSR B0656+14 in the / and 
R bands basically completes a comprehensive broadband pho- 
tometric study of both pulsars in the NIR through FUV range 
started about a dozen years ago. As a result, the two pulsars 
have been detected at a high significance level practically in 
all standard optical bands generally used for the photometry 
of ordinary stars to obtain the spectral energy distributions of 
their radiation in a wide spectral domain. Except for the most 
bright and younger Crab and Vela pulsars, Geminga and PSR 
B0656+14 stand up to be the most studied up to date radio pul- 
sars in the optical range. 

4.1. Spectral energy distribution in the optical range 
4.1 .1 . Extinction corrections 

Interstellar color excesses E{B - V) along the lines of sight to 
Geminga and PSR B0656-(-14 have not yet been investigated 
accurately. This requires additional careful optical studies of 
suitable stars from the pulsar neighborhoods. So far there are 
only indirect estimations of the extinction based on observa- 
tions in other spectral domains (e.g., Pavlov et al. .1996 Kurt 
et al. [Tp98 Kurt et al. 2001). Closeness of both pulsars 
to the Galaxy plane (see Table ^ suggests that their optical 
fluxes are to be affected by the foreground extinction from the 
Milky Way. On the other hand, the reddening cannot be too 
strong due to their proximity to the Earth. This is confirmed 
by relatively low DMs and hydrogen column densities Nh ob- 
tained from the radio and X-ray observations. The conserva- 
tive Nh values are a; 1.1 x lO^'^cm"^ for Geminga based on 
the ROSAT and ASCA X-ray spectral observations (Halpern 
& Wang [19971 Jackson et al. l2U02b and =sl.4 x lO^^'cm ^ for 
PSR B0656-H 14 based on the ROSAT, ASCA and EUVE obser- 
vations (K2001). Using these values and an empirical relation 
Nh/E(B - y) = 4.8 X lO^^cm-^mag-i applicable for the Milky 
Way (Bohlin et al. \WM one can obtain E(B - V) of 0.023 
and 0.03 for Geminga and PSR B0656-I-14, respectively. These 
values are consistent with those suggested by previous studies. 
Apparently smaller value of E{B - V) for Geminga is compati- 
ble with the fact that it is by about 100 pc closer to us than PSR 
B0656+14. We applied these excesses to calculate dereddened 
fluxes in a standard way making use of mean Galactic extinc- 
tion curve (Cardelli et al. 1989 1 and the de-reddening param- 
eter R{V) = 3.1. The results are presented in Tables |5] and |6l 
For both pulsars the dereddening corrections are within 1 cr of 
the flux uncertainties in the optical-NIR range, while they are 
significant in the UV. 

We note that spectral analysis of recent Chandra and XMM 
X-ray observations with previously used standard spectral 
models give systematically higher Nh values than those of the 
ASCA and ROSAT: *2.4xl02°cm"2 ^j. Geminga (Kargaltsev 
et al. 12005 J and «;4.3xl0-°cm-2 for PSR B0656+14 (De Luca 
et al. l2005i . The latter one corresponds to E(B - V) ~ 0.09 
which is comparable with the extinction throughout the entire 
Galaxy in this direction (Schlegel et al. .1998} and hardly can 
be relevant for the nearby pulsar. Nh for Geminga appears to be 



also overestimated and we skip out these implausibly high val- 
ues '° from our consideration and use the above ones provided 
by early ASCA, EUVE, and ROSAT data analyses. We note, 
that Nh was also fixed at the ROSAT value in the recent analy- 
sis of the Geminga X-ray spectrum based on the most complete 
data set taken with the XMM (Caraveo et al. 2004). 

4.1 .2. Broadband spectra from NIR tlirougli UV 

The dereddened broadband spectral energy distributions of the 
two pulsars are shown in Fig.0 Despite a factor of 2-3 differ- 
ence between the mean flux levels of the two objects and higher 
flux uncertainties in case of Geminga, the spectral shapes of 
the optical emission of both middle-aged pulsars are remark- 
ably similar to each other Both are clearly nonmonotonic with 
apparent dips in the F430W and FllOW bands and excesses 
centered in the V-I range. PSR B0656-(-14 has also a signifi- 
cant flux increase towards lower frequencies in the NIR range. 
Similar increase is seen in the spectrum of Geminga but it is 
less significant due to the absence of observations in the K 
band (or in its HST analog F187W). There is also flux in- 
crease towards higher frequencies in the UV range which is 
most prominent for Geminga where the FUV data are available. 
The spectroscopic data recently obtained for PSR B0656+14 
in the NUV (Shibanov et al. 2005 ) and for Geminga in the 
FUV (Kargaltsev et al. l2005t are also shown in Fig. 7 for com- 
parison. They are compatible with the broadband data, but the 
flux uncertainties (not shown) in spectral bins are considerably 
larger than photometric ones, and apparent absorption/emission 
features in the spectra have a low significance of ^ 2cr. 

Resulting spectral energy distributions cannot be fitted by 
a single power-law, as it is possible for the young Crab pul- 
sar in the whole observed range (Sollerman '2003" SoUerman 
et al. i2000:) . This suggests a multicomponent structure of the 
emission of the two middle-aged pulsars and may reveal spec- 
tral evolution of the optical emission with pulsar age (K2 001I 
Shibanov et al. .2003. Zharikov et al. .2004J . Some of the com- 
ponents can be identified from the comparison with the X-ray 
data. 

4.1.3. Tliermal and nontliermal spectral components 

For PSR B0656+14 the flux increase with frequency in the UV 
can be explained by a contribution of the Rayleigh- Jeans part 
of the thermal blackbody-like spectral component (BB) from 
the whole surface of the NS strongly dominating the pulsar 
emission in soft X-rays. The flux increase towards the infrared 
appears to be consistent with the low energy extension of the 
power law component (PL) seen in the high energy tail of the 
X-ray spectrum. The unabsorbed PL and BB components and 
their sum are shown by dashed lines in Fig. 7. Their parame- 

The overestimation of Nh in both cases is likely due to the poorly 
calibrated instrument responses in the energy range of S 0.3 keV 
(e.g., Kargaltsev et al. 2005 i which is most crucial for Nh estimates. 
This also leads to lower temperatures of thermal emission components 
and to higher spectral indices of nonthermal components derived from 
spectral fits. 
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Fig. 7. Dereddened NIR-UV broadband spectra of PSR B0656+14 (left) and Geminga (right) summarizing the best up to date 
quality observations obtained with different telescopes, instruments and filters, as indicated in the plots and listed in Tables 5 and 
6. Note a nonmonotonic structure and remarkable similarity of the broadband spectra for both pulsars. Thin solid curves show 
the spectroscopic data obtained for PSR B0656+14 in the NUV (Shibanov et al. 2005 1 and for Geminga in the FUV (Kargaltsev 
et al. I2005> . Dashed lines are the low energy extensions of the unabsorbed soft blackbody (BB) and power law (PL) spectral 
components derived from the combined BB-i-BB+PL spectral fits in X-rays. Parameters of the fits are taken from Koptsevich et 
al. (2001 1 for PSR B0656+14 and from Kargaltsev et al. (2005 1 for Geminga. The second, hard BB components are negligible 
for both pulsars in the NIR-UV range and not shown in the plots. Dotted and dot-dashed lines in the right panel show the PL fit of 
the Geminga optical data in the range from F160W through F430W bands and its sum with the BB, respectively (see Sect. 4.1.3 
for details). 



ters are taken from H2001I the BB temperature Tbb~OM MK 
and the radius of the emitting area /?bb^7.8 km at the dis- 
tance d~2^S pc; the PL spectral index a'vSi-0.45 and normal- 
izing constant at 1 keV is r;2.8x10"^ ph cm"^ s ' keV '; 
Nh~ 1 .42x 10^° cm-2 ^jjg ^^^^ E(B-V) in the plot are 

consistent with each other as discussed above. Uncertainties of 
the X-ray spectral fit extension are shown by a dashed error- 
bar As seen from the left panel of Fig. 7, the BB-hPL extension 
qualitatively fits the spectrum of PSR B0656H-14, albeit the PL 
component dominates in the most part of the observed optical 
range. This is also confirmed by detailed analyses of the pul- 
sar pulse profile and the spectrum in the NUV, where the BB 
contribution was estimated to be ^ 30% (Shibanov et al. 120051 . 

For Geminga the situation is more complex. The param- 
eters of the X-ray extensions in this case are taken from 
Kargaltsev et al. (|20^ : Tbb~0.52 MK, Rbb~9.5 km at 
d^2QQ pc; the PL index ay~-Q.56 and normalizing constant at 
1 keV is «i5.5xlO"^ ph cm"^ s"' keV"'; Nh was fixed during 
the X-ray spectral fit at the value of 1.1x10^*^ cm"^. As seen 
from the right panel of Fig. 7, the PL and, therefore, the sum 
BBh-PL (it is outside the upper plot frame) overshoot the opti- 
cal fluxes of Geminga by a factor of 100. This fact has been 
noticed by Halpern and Wang ( 1 997j and confirmed later on by 
Jackson et al. | |2002> . Komarova et al. J2003> . and Kargaltsev 
et al. J2005> . At the same time, the BB extension appears to 
fit qualitatively well a steep flux increase of Geminga towards 
higher frequence range. 



Various possibilities to resolve the discrepance between the 
low energy extension of the Geminga X-ray spectrum and the 
observed optical fluxes have been discussed by Kargaltsev et 
al. ( 2005 1, including possible systematic errors due to instru- 
ment response uncertainties, large statistical uncertainties of 
the X-ray PL component, its contamination by the hard BB 
component from the pulsar polar caps, and possible existence 
of a spectral break in the spectrum of the nonthermal compo- 
nent somewhere between the optical and X-ray ranges. 

The similarity of the optical spectral shape of Geminga with 
that of PSR B0656-I-14 suggests that its low frequence tail is of 
nonthermal origin and can be fitted by a power law. Such a fit 
performed for the F160W-F430W range is shown by a dotted 
line in the right panel of Fig. 7. Its sum with the BB extension 
(dot-dashed line) qualitatively fits the whole set of the data in 
the NIR-UV range, as it is in case of the BB+PL extension 
for PSR B0656-I-14. The slope of the optical PL, q'v~-0.26, 
is almost twice smaller than that in X-rays, ayX-Q.56, favor- 
ing the interpretation with the spectral break in the nonther- 
mal optical-X-ray emission of Geminga. Changing the spectral 
slope of the nonthermal component from almost a flat in the 
optical to a significantly negative in X-rays is not unusual and 
has been observed in the emission of the young Crab pulsar 
(Kuiper et al. 2001| SoUerman .2003i) . older Vela (Shibanov et 
al. 2003 i Romani et al. l2005> . and hkely in two old pulsars PSR 
B1929-H10 and PSR B0950+08 detected in the optical (Pavlov 
et al. ,.1996. Mignani et al. 2002 ; Zharikov et al. .2002., 2004 j and 
X-rays (Becker et al. .2004 and 2005. Zavlin & Pavlov l2004bt . 
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The nonthermal emission of very old isolated millisecond pul- 
sars studied in the optical and X-ray ranges shows even a 
much stronger slope break between these ranges (Koptsevich 
et al. 2003, Mignani & Becker 2004). We discuss the possible 
position of the break for Geminga below in Sect. 4.2. 

4.1 .4. Third optical spectral component? 

A wide F555W-R excess in the emission of PSR B0656-i-14 
cannot be explained by only the combination of the Rayleigh- 
Jeans and power-law components. It exceeds by ~ (3 - 5)cr 
a power-law "continuum level" supposed to be seen in the 
F43QW and / bands. The excess may reveal a new, not yet 
identified, third spectral component in the emission of PSR 
B0656-H14. 

The nature of the excess is unclear. One can speculate that 
it is produced by unresolved emission lines from a compact 
PWN around the pulsar interacting with interstellar environ- 
ment polluted by heavy elements from the SN explosion asso- 
ciated with the pulsar. Strong Olll (5007 A), Sll (67 17 A) and 
other fainter nebular lines, whose wavelengths are within the 
spectral range of the excess, are typical for the spectra of bright 
and extended PWNe formed around young pulsars like Crab. 
In case of PSR B0540-69 a strong background emission in Olll 
(5007 A) line from the SNR core contaminates the pulsar pho- 
tometric flux in the F555W band and results in a false flux 
excess exceeding by an order of a magnitude the real flux from 
the pulsar (Serafimovich et al. 12004 j . A distant outer shell of 
the PSR B0656+14 SNR has been detected in X-rays (Nousek 
et al. ll981l Thorsett et al. 2003 1. However, nothing is known 
on the presence of the nebular lines in the pulsar vicinity, al- 
though they could exist if the pulsar PWN marginally detected 
in the optical (Fig. 6) and X-rays (Marshall & Schulz l2002l l is 
real. 

Alternatively, the excess can be an internal property of the 
pulsar emission. It can be produced by ion or electron cyclotron 
features formed in magnetospheric plasma of a strongly mag- 
netized NS. The features are expected to be broad enough due 
to the magnetic field inhomogeneity at the magnetospheric spa- 
tial scale. Current data do not allow us to discriminate between 
the two possibilities. To do that a deep narrow-band imaging 
of the pulsar field with a standard set of narrow-band filters 
centered at nebular lines and/or spectroscopy of the pulsar are 
necessary. 

Geminga is fainter and its fluxes are determined with larger 
errors. The apparent I-R excess in its spectrum is much less sig- 
nificant than for PSR B0656-I-14. However, the similarity of the 
spectral shapes of both pulsars provides an additional argument 
that this excess is real. This is also supported by the presence of 
the X-ray nebulosity around Geminga and its possible optical 
counterpart (Fig. 5), which may be resposible for the excess. In 
this case, one can suggest the same two possibilities to explain 
its nature as has been discussed above for PSR B0656-I-14. 

We note, that with the new optical and NIR data the optical 
spectrum of Geminga cannot be described by a two component 
model (Jacchia et al. 1999.) combining only an ion cyclotron 
emission line from a hot plasma near the polar caps of the pul- 



sar and the thermal radiation from the cooling surface of the 
NS. The nonthermal PL spectral component has to be included 
into the interpretation and it strongly dominates the emission at 
longer wavelengths. 

4.2. Multiwavelength spectra 

In Figure|S]we show unabsorbed rotation-phase integrated mul- 
tiwavelength spectra of both pulsars using available data from 
the radio through gamma-rays. The radio data for Geminga 
are kindly provided by V. Malofeev (private communication). 
There are some doubts on the radio detection (Burderi et 
al 1999) and the data are represented mainly by upper lim- 
its. For PSR B0656+14 the radio data are from i dnoTI The 
optical data are from Tables 5 and 6. The parameters of the X- 
ray spectral components are the same as in Fig. 7. For both 
NSs we add here also the hard BB components believed to 
be from the pulsar polar caps. They contribute significantly 
to the resulting spectra only near ~1 keV where the soft BB 
components meet the PL tails. Parameters of the hard BBs are 
from K'MJTfor PSR B0656-H14 (R^O.52 km at d~288 pc, 
ra;1.65 MK) and from Kargaltsev et al. (20P5) for Geminga 
(7;^0.046 km at d~200 pc, Tk2.32 MK). AU X-ray compo- 
nents are extended to the optical. The gamma-ray data are 
from Ramanamurthy et al. jl996> for PSR B0656-I-14 and from 
Jackson et al. (2002 RXTE), Macomb and Gehrels ) fT^ 
OSSE, COPMTEL, E GRET) ", Akerlof et al. ) fT993l Whipple) 
and Vishwanath et al. ( 119931 Tata group) for Geminga. 

The multiwavelength spectra underline the fact that a 
smooth connection of the optical and X-ray nonthermal and 
thermal spectral components observed for PSR B0656-)-14 
is problematic for Geminga. Geminga is obviously under- 
luminous in the optical range as compared with PSR B0656-)-14 
and with the low-energy extension of its power-law component 
from X-rays. This suggests changing the spectral slope of the 
nonthermal component of Geminga from a significantly neg- 
ative in X-rays (q'v~-0.56, dashed line) to almost a flat (oy^- 
0.26, dotted line) in the NIR-optical range. A simple extrapola- 
tion of the NIR-optical PL fit to X-rays shows that the spectral 
break in the nonthermal optical-X-ray emission may occur near 
1 keV where the hard BB component dominates the total flux 
(bottom panel of Fig. 8). 

A common origin of the optical and X-ray nonthermal com- 
ponents independently on pulsar age and the presence of spec- 
tral breaks was proposed by Zharikov et al. ( 2004"; '2005') based 
on a strong correlation of the optical and X-ray nonthermal 
luminosities for a set of ordinary pulsars detected in the op- 
tical and X-ray domains (see also Zavlin & Pavlov '2004b1. 
Analysis of the phase aligned optical and X-ray pulse profiles 
would be useful to verify this hypothesis for Geminga, as it has 
been partly done for PSR B0656-I-14 (Shibanov et al.|^05j. 
Detection of the optical pulsations of Geminga at a higher 
S /N than it is provided by a tentative detection by Shearer et 
al. (I1998> is necessary for this goal. 

The multiwavelength nonthermal spectrum of Geminga 
suggests also at least three additional spectral breaks in the 

' ' ^jssc.gsfc. nasa.gov/docs/cgro/gamcat/ 
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Fig. 8. Unabsorbed multiwavelength spectra of the two pulsars using available data from radio through gamma-rays, as notified 
in the plot. X-ray spectral components are extended to the optical range, while the optical one - to X-rays. 



high energy range: in hard X-rays, in gamma-rays, and in the 
TeV range. They form a "knee-like" spectral structure of the 
high energy emission. Similar "knee" is seen in the high en- 
ergy spectrum of the Vela-pulsar (e.g., Shibanov et al. 120031 



Romani et al. 12005 1 and references therein). The optical part of 
the Vela spectrum suggests also a "degenerate knee" between 
the optical and X-rays while its NIR and X-ray data from the 
nonthermal spectral tail can be fitted by a single power-law 
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(Shibanov et al. 120031 ). It is not clear whether the optical-X- 
ray and X-ray-gamma-ray knees are correlated or physically re- 
lated to each other. The Crab pulsar does not have such a strong 
spectral breaks and knees in the high energy range, while it has 
a strong spectral break between the optical and X-rays (Kuiper 
et al. 1200 1 Sollerman 2003 1. The same appears to be true for 
PSR B0540-69 (Serafimovich et al. 2004). Various slopes and 
numerous breaks likely suggest a multicomponent structure of 
the emission generated by different mechanisms or by different 
populations of the relativistic particle, or both. However, multi- 
wavelength data for pulsars are too scarce to make any definite 
conclusion. For instance, PSR B0656-H14 appears to have no 
significant spectral break between the optical and X-rays, while 
almost nothing is known on its high energy spectrum except it 
is likely softer than that of Geminga. 

The problem of the spectral difference of the pulsars out- 
side the radio range is far from its final solution. A question, 
whether this difference simply caused by different pulsar emis- 
sion beaming and geometry factors or by significantly differ- 
ent physical conditions in the pulsar magnetospheres, is not 
answered yet. Studying the pulse profiles in different spectral 
ranges and phase resolved spectroscopy can help to solve this 
problem and, particularly, to understand whether the optical 
and X-ray nonthermal photons are generated at the same site 
and by the same population of relativistic particles accelerated 
in the pulsar magnetospheres. If the latter is true, then one has 
to understand why the distribution function of the particles in 
some cases reveals a break at an intermediate energy, which is 
reflected in the observed spectra. 



4.3. Possible optical PWNe 

The sizes of axisymmetrical equatorial torus-like (plus polar 
jets) PWNe around the young Crab, PSR B0540-69, and Vela 
pulsars are roughly scaled as E^-^, where E is the spin-down 
luminosity of a pulsar (Helfand et al. 2001 Pavlov et al. 2001 
Serafimovich et al. 120041 Weisskopf et al. .2004j . This is in 
agreement with the simplest model of a termination wind shock 
forming around subsonic-ally moving energetic pulsars due to 
the interaction of their relativistic winds with the ambient mat- 
ter (Kennel & Coroniti 1984 1. The PWNe of the Crab and its 
twin PSR B0540-69 are extended approximately to about 1 pc 
from the pulsars and both have similar sizes in the optical and 
X-rays (e.g., Serafimovich et al. 2004). Using this value and 
the above scaling law one can estimate that PWNe around the 
middle-aged Geminga and PSR B0656-I-14 have to be about 
a factor of 100 smaller than in the Crab case, i.e, a few arc- 
seconds at the distances to the pulsars of about 200-300 pc, 
provided the external density and pressure conditions are not 
much different from that of Crab and are appropriate to form 
the subsonic torus but not a supersonic bow shock. The sizes 
of the faint extended structures around both middle-aged pul- 
sars marginally detected in our deepest Subaru images are in a 
reasonable agreement with these expectations. This is a start- 
ing point to discuss these structures as possible candidates for 
PWNe. 



So far PWNe of a different nature and shapes have been 
detected in the radio, optical and X-rays around a dozen of iso- 
lated NSs (e.g., Gaensler 2005'). Young NSs typically show the 
torus/jet-like PWNe, while the PWNe around old NSs com- 
monly have a bow shock structure. Changing of the shock type 
with age is supposed to be due to a transformation of exter- 
nal conditions from the subsonic to supersonic ones when the 
density and the pressure within the expanding host SNR de- 
crease to values of a typical ISM. A characteristic age of a host 
SNR when it dissolves in the ISM background roughly coin- 
cides with a middle-age epoch of the NS evolution, i.e., 10^ - 
10'' yr. 

Following this scenario one can expect to see an inter- 
mediate PWN morphology, i.e., a mixture of torus and bow- 
shock, around middle-aged pulsars. However, until now there 
are only two middle-aged pulsars with a faint, few arcsec- 
onds PWNe apparently detected only in the radio around PSR 
B0906-49 (Gaensler et al. [T9981 and only in X-rays around 
PSR J0538+2817 (Romani & Ngl2003l. Both cases need ad- 
ditional observations in different wave bands to confirm the na- 
ture of the extended objects around the pulsars and to study 
their morphology. In this respect, adding to this poor sample 
any detection of PWNe, even a tentative, would be valuable. 

In contrast to a bright Crab PWN, whose contribution to the 
total pulsar -i-nebula luminosity is about 95.5% and 99.8% in 
X-rays and in the optical, respectively, the PWN-like structures 
around our middle-aged pulsars are faint and may contribute 
only a few percents to the total emission budget. However, 
this is compatible with the observational fact that the torus 
PWNe loose their brightness with age much sooner than the 
pulsars. For instance, the PWN contribution for an intermedi- 
ate age Vela-pulsar (10'* yr) in X-rays is about 75% (Pavlov 
et al. 2001). It is likely to be much smaller in the optical 
range where the Vela-PWN has not yet been clearly identi- 
fied (Shibanov et al. 2003; Mignani et al. 2003). A faint X- 
ray PWN tentatively detected around a middle-aged pulsar PSR 
J0538-H2817 (6 x 10^ yr) provides only about 2% of the pulsar 
flux in X-rays (Romani & Ng 2003 ). This is in agreement with 
what is seen in our case of Geminga and PSR B0656+14. 

Another factors suggesting that the extended optical struc- 
tures around our pulsars are likely to be PWNe are their mor- 
phology and orientation to the vectors of the pulsar proper mo- 
tion. The planes of the equatorial torus structures of young 
pulsars are typically arbitrary inclined to the line of sight and 
seen as elliptical extended objects. For the Crab, Vela, and, 
possibly, PSR B0540-69 pulsars the torus symmetry (pulsar 
spin) axis and PWN polar jets are aligned with the vector of 
proper motion of the pulsar, which may constrain the nature 
of the kick at pulsar birth (Serafimovich et al. 20CP5. Both 
properties are likely seen for Geminga and PSR B0656H-14 
where the extended optical structures (possible inclined tori) 
are aligned perpendicular to the pulsar proper motion, although 
any jet/counter-jet structures are not visible in our images. 

Finally, possible pulsar nebulae around Geminga and PSR 
B0656+14 have been detected in X-rays, albeit at different sig- 
nificance levels. The sizes of a less significant circle X-ray neb- 
ulosity around PSR B0656-H14 (Marshafl & Schulz lMHt are 
consistent with what is seen in the Subaru image. Geminga 
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shows apparently more significant long (up to ~2') X-ray 
tails well aligned with the source proper motion (Caraveo 
et al. I2003I I. The tails evidence that the PWN of the older 
Geminga has transformed from the subsonic torus to the su- 
personic bow shock stage. This is also supported by a twice 
higher transverse velocity of Geminga, by absence of its host 
SNR, and, hence, likely smaller density and external pressure 
of the ambient matter than in case of PSR B0656+14. The lat- 
ter one is younger and still sits inside its host SNR (Thorsett et 
al.p003l. Note also, that its DM is significantly and Nh is ap- 
parently higher than in the Geminga case, favoring a denser en- 
vironment for PSR B0656+14. Consequently, an arc-like shape 
of the faint extended structure at the Geminga optical image 
suggests that it is likely to be a counterpart of the head of the 
bow-shock (but not a bright part of a PWN torus) whose tails 
are detected in X-rays. 

Thus, the morphology of both optical nebulosities tenta- 
tively detected around the two middle-aged pulsars are roughly 
compatible with what is expected from their parameters and pa- 
rameters of the environments, as well as with what is seen for 
each pulsar in X-rays. Confirmation of their reality at higher 
significance level and with better spatial resolution in the opti- 
cal and X-rays would be valuable for understanding the physics 
of PWNe and their evolution with age. 
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